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ABSTRACT 

An effective data interface between a Victoreen PIP-400 
pulse-height analyzer and a Hewlett-Packard 9810A calculator 
was designed, built, and tested. A calculator program was 
written which enabled a research group studying laser 
propagation in the marine boundary layer to conduct rapid, 
local processing of scintillation and extinction data. 
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I. 



INTRODUCTION 



The potential utility of laser technology in the areas 
of line-of-sight communication, target tracking, and wea- 
ponry makes an understanding of the propagation of laser 
beams near the air-ocean interface imperative. The effects 
of a turbulent atmosphere on laser transmission were theor- 
etically modelled by Tatarski [Ref. 1] et. al. and, subse- 
quently, extensive field research in the terrestrial boundary 
layer has substantiated the bulk of the theory. Research 
data in the marine boundary layer has been scarce and 
inconclus ive . 

Atmospheric effects on laser propagation may be divided 
between laser- induced phenomena and those caused by meteoro- 
logical characteristics independent of the laser's presence. 
The first classification applies only to lasers of sufficient 
power to heat the atmosphere (thermal blooming, upwind 
drift) . The second is a universal grouping affecting the 
whole electromagnetic spectrum regardless of power. The 
atmospheric effects include beam refraction, beam spread 
(or degradation of resolution) , absorption, and scattering. 
Refraction may be further categorized according to the 
dimensions of the atmospheric whirls. Large-scale vortices 
cause path deviation (beam wander) while small-scale vortices 
cause time- intensity variations (scintillation) and degrada- 
tion of resolution. 
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The Naval Postgraduate School, Monterey, has been involved 
over the last four years in a continuing effort to establish 
the data base necessary to confirm or modify existing theory 
as applied to the air-ocean interface. Two research groups 
have operated simultaneously but independently; the optical 
research group and the atmospheric turbulence group have 
attempted to correlate determinations of (the refractive 
index structure constant) through the respective analysis 
of light intensity distributions and meteorological data. 

The optical research group has concurrently investigated 
scintillation, extinction (scattering and absorption) and 
the atmospheric modulation transfer function (a measure 
of beam spread and wander) . Theses on the effects of aperture 
averaging, determination of an atmospheric extinction con- 
stant, and determination of an atmospheric MTF for a black- 
body source were completed over the same period as the 
author's [Refs. 10, 11, 12]. 

A major problem for the optical group was the inordinate 
amount of time necessary to process scintillation and extinc- 
tion data. The methods employed depended heavily on operator 
analysis of raw data or extensive manual input to computers. 
While the processing theory being applied was appropriate, 
the time penalty and the inherent high probability of error 
had to be minimized. The author concentrated his efforts 
in this area employing only the basic results of laser 
propagation theory . 
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II. PROJECT BACKGROUND 



Various members of the Naval Postgraduate School's 
faculty have been conducting research on atmospheric pro- 
pagation phenomena since 1972. The first experiments inves- 
tigated the atmospheric MTF and scintillation in the basement 
corridor of Spanagel Hall. In the scintillation experiment 
a laser beam reflected from a flat mirror placed at the 
opposite end of the 137 meter passageway and entered a pin 
diode detector located near the laser source. The detector 
output, proportional to the incident beam intensity, was 
passed through a log- amplifier and recorded on magnetic tape. 
The data was later periodically sampled, digitized, and 
accumulated by a pulse-height analyzer; the resultant dis- 
tribution was Gaussian, as theory had predicted. Refractive 
index structure constants were calculated from photographs 
of the distributions by measuring the curves ' full width at 
half maximum height. The MTF measurement used a six inch 
parabolic mirror mounted as a Newtonian telescope. The 
image from a laser at the other end of the passageway was 
scanned across a slit using a cam-driven rocking diagonal 
mirror. The output of a photocell was fed to a boxcar inte- 
grator to generate time-averaged data, which was then recorded. 
Turbulence in the hall was somewhat controllable and detec- 
table reduction of MTF with turbulence was observed. 
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In late 1972 a group of interested researchers discussed 
possible environmental experiments. The first such experi- 
ment was completed 10 January 1973 by faculty, staff, and 
students and involved shining a low-power (mW) red laser 

O 

(6328 A) across a 4.3 kilometer range at the southern end 
of Monterey Bay. The research vessel Acania was anchored 
at the approximate midpoint of the beam path to take atmos- 
pheric turbulence measurements with onboard sensors. In a 
second experiment at San Nicholas Island (3 March 1973) low 
power laser beams were transmitted from the northwest tip 
of the island to the Acania anchored offshore. The qualified 
success of these experiments completed with equipment already 
owned by the school was encouraging enough to prompt the 
submission of a research proposal to the Naval Ordnance 
Laboratory (White Oak, Md.) for fiscal year 1974. 

Experimental apparatus has changed significantly since 
the Navy lab started sponsoring the research. A Cassegrainian 
telescope system with an eighteen inch parabolic mirror and 
an elastically mounted scanning mirror (driven electromag- 
netically by a shaped waveform) was used in the most recent 
MTF measurements. An enclosed trailer was built to transport 
and protect the telescope and associated equipment; aluminum 
support legs that extended to rest directly on the ground 
negated the previously persistent problem of telescope motion 
due to operator movement. More recently a bus has been used 
which holds most of the on-shore optical and meteorological 
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equipment and employs the same independent telescope suspen- 
sion idea. Ship motion originally made it difficult to keep 
even a diverged laser beam on the telescope long enough to 
collect maningful MTF data. A gyro-stabilized tracking system 
was designed to go aboard the Acania. At first it was manned 
by two operators who sight-aimed a gyro-mounted laser with 
hand-operated servo joy sticks; one controlled horizontal 
and the other the vertical train. Little useable data was 
obtained with this system because of problems with human 
coordination and fatigue. The solution was a quadrant detec- 
tor servo system which will lock on a laser beam once the 
operator has aligned it with the target. Shipboard scin- 
tillation detectors are also mounted on the gyro frame and 
thereby maintain a constant relative position in the incident 
beam path. The gyro has been housed in an aluminum shed, 
when aboard the Acania, to protect both equipment and operator 
from the elements . 

Scintillation data was best taken aboard the Acania; a 
detector with a wide angle of acceptance could continuously 
sample a diverged beam from ashore more readily than the 
reverse because of the effect ship motion would have on the 
direction of propagation of a beam originating on board. 

The single-detector system originally employed required 
excessive attention during data runs because light from 
environmental sources refracts into the beam path and can ' t 
be distinguished from that of the monitored laser. Since all 
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The bus with 
telescope in 
position. Transit 
to the right was 
used in range-to- 
ship calculation. 

In background stands 
meteorological tower 
supporting sensors. 




Naval Postgraduate 
Research Vessel Acania 
with meteorology gear 
deployed. The 
aluminum shed which 
shelters the gyro 
is aft of the 
pilot house. 
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Right side of 
telescope showing 
fan^-scan unit, 
laser, detectors, 
tracking unit, 
and associated 
electronic gear. 



Gyro-stabilized 
tracking unit. The 
detectors are 
mounted on top and 
the lasers are 
hidden from view 
below. 
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Oscilloscopes , 

Nuclear Data Analyzer, 
log-amplifier, and 
demodulator; here 
shown set up for 
meteorological data 
processing on the 
bus . 



Dry lab on the 
Acania with 
scintillation/ 
extinction data 
processing gear: 
PIP-400, PAR 
amplifiers , 
demodulator, log- 
converter, HP9810A 
calculator and 
plotter. Oscillo- 
scope & freq 
analysis tape 
recorder . 
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photocell data passed through a log-amplif ier a skewed 
distribution would result. In an attempt to negate this 
effect the beam had to be interrupted ten seconds out of 
every minute to zero-out detector voltage generated by back- 
ground radiation. A dual-detector system was built with the 
idea that the output of one pin diode detecting only background 
could be subtracted continuously from that of another unit 
aimed at the laser beam. 

Since the detectors (supposedly identical) produced 
different proportional outputs extensive voltage biasing 
was necessary. An even greater problem was that ship motion 
continuously changed the position of the detector relative 
to the beam center. Subsequently another single-detector 
system was designed, this one employing a chopped laser source 
and a demodulator. The chopping was accomplished at the source 
by reflecting the laser light first off of a cylindrical 
mirror which fanned the beam horizontally and then off of 
a flat mirror vibrating at approximately two kilohertz which 
made the beam scan vertically. The demodulator unit took 
the detector output and electronically subtracted the minima 
(background) from the maxima (laser & background) and sent 
the difference to the log-converter. Detector sensitivity 
was recently increased by replacing the pin diode with an 
avalanche detector. An FET follower converts the high- 
impedance output of the detector to a low impedance which 
can be transmitted a greater distance through coaxial cable. 
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This has allowed all data processing gear on the Acania to 
be positioned below-decks in the dry lab. 

Processing of scintillation data was initially best 
accomplished through physical measurement of photographs 
of resultant distributions displayed on the PIP-400 pulse 
height analyzer. Attempts to record and subsequently analyze 
data with integral methods on a computer suffered from poor 
digitizing. Digital output could be obtained directly from 
the PIP-400 by a fast printer or from a Nuclear Data 128- 
channel analyzer by teletype. An attempt was made to fit a 
Gaussian curve to the data using least-squares method. This 
analysis proved highly successful when extraneous voltages 
did not skew the experimental distribution but was painfully 
slow because data had to be manually input to the computer 
point by point. The most recent improvement has been the 
direct interfacing of the PIP-400 with a compact computer 
programmed to process scintillation and extinction data on- 
scene within six minutes of sampling. This was completed 
by the author and is the subject of this thesis. 
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III. THEORY 



A. SCINTILLATION 

The theory of optical propagation in a turbulent atmos- 
phere has been extensively documented and those portions 
employed by the Naval. Postgraduate School optical propaga- 
tion group have not changed significantly since the project 
began. A thesis submitted by B. C. Haagensen [Ref. 3] in 
June, 1973 contained an excellent outline of the pertinent 
theoretical derivations and an extensive bibliography. The 
author does not intend to list or verify familiar theoretical 
developments but rather to present the significant end 
products vital to data processing. 

Theory states that atmospheric turbulence is a random, 
non-linear, three-dimensional phenomenon. A turbule is 
ideally assumed to be an isotropic, homogeneous pocket or 
swirl of air with a particular index of refraction. The 
inertial subrange delineates a region of turbule size in 
which refraction is theoretically described; it is bounded 
on one side by turbules so small (a few millimeters) that 
viscous dissipation of energy causes their immediate collapse 
and on the other by eddies so large (order of meters) that 
they are no longer homogeneous and energy is added to the 
turbulence spectrum. In the inertial subrange large turbules 
form and transfer energy to successively smaller turbules 
with minimal losses through viscous dissipation. 
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Tatarski [Ref. 1] described the three-dimensional power 



spectral density of turbulence, $ n (K), with the equation: 



obviously parameterizes the total amount of energy 
in the turbulence. The random distribution of turbulence 
cannot be forecast and must be described statistically. A 
beam is modulated in a multiplicative manner when it passes 
through successive turbules of differing index of refraction. 
It follows that the natural log of the beam intensity is 
altered additively and should follow a normal probability 
distribution as a result of the central limit theorem. C XT 
is directly related to the structure function of the refrac- 
tive index (the mean square difference of the index of 
refraction at two points a distance d apart) by: 



* (K) = 0.033 C. 
n 



’N 



2 K ' 11 / 3 




where : 



G., = refractive index structure constant 



2 7T 

K = spatial wave number = — 
t = coherence length of a turbule 



2 

K = = lower limit of inertial subrange 

m l n 0 



0 



D n (d) = < [n (r) - n(r+d)] 2 > = C N 2 d 2// 3 



20 



After lengthy derivation and a number of assumptions [Ref. 
4, pg. 10], can be seen to be directly dependent on the 
variance of the log intensity distribution by: 



where 



C N = 1.42 K- ?/12 Z- 11/12 c tl 



2 2 

= <[ln(I) - ln<I>] > = variance 



= optical propagation path length 



can also be calculated by a statistical determination 

2 

of the temperature structure constant, C T , since tempera- 
ture is the dominant factor in a meteorological formulation 
of the index of refraction: 



C N = (79.0 ^ x 10 6 ) C, 



where 



P = mean pressure in millibars 

T = mean temperature in K 

„ 2 <[T(r) - T (r+d) ] 2 > 

C T = d 2/3 
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A major objective of the project is to study the 
correlation between C N values obtained optically and 
meteorologically . 

B. GAUSSIAN MODELLING 

According to theory, periodic sampling of a laser beam 
propagating in the inertial subrange should approximate a 
lognormal distribution; if the natural log of the detector 
voltage (proportional to beam intensity) defines abscissa 
values and the ordinate marks the number of samples returning 
a specific voltage the resultant plot should be of the form: 




2 



where : 



C is a constant 



x is natural log of detector voltage 



x 



is the average natural log of detector voltage 



2 2 
a is the variance (C^ 



Taking the natural log of each side: 



2 



2 ^ 

ax + bx + c 



In y = - 



x 



2a 



2 + *7 + (lnC - 7T } 

a 2a 
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where : 



1 



a 



2 



2a 



b 



x 



2 



a 



c = In C 




2a 



2 



Experimentally the natural log of detector voltage axis 
is divided by the PIP-400 into two hundred channels of equal 
width and the y axis is the number of counts (or times) the 
equivalent detector voltage (or beam intensity) was detected 
over the sampling period. The number of counts in each chan- 
nel after a data run can be assumed to be Poisson-distributed; 
for one measurement the 68% confidence interval is y ± /y. 
Since the natural log of y is taken to obtain the quadratic 
form, the uncertainty in ln(y) becomes; 



So, for the log distribution, the channels with the largest 
numbers of counts have the narrowest confidence intervals 
(smallest uncertainty) and should be weighted more heavily. 

The effect of narrowing confidence intervals with increasing 
counts is easily illustrated by plotting them on log paper. 

In the original program it had been decided to weight the data 
inversely with the uncertainty, i.e. multiply by /y. 



d [In (y ) ] 



dy _ /y _ 1 

y y jz 



/y 
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THE LOGNDRMRL 
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The ln(y) distribution resulting experimentally can be 
approximated with a polynomial by means of least-squares 
fitting. The method seeks to minimize 

199 

Z w . [ In ( y . ) - ax. 2 - bx. - c] 2 
i=0 1 1 1 1 

where: 



x^ = channel number i 
y^ = counts in channel i 

w^ = weighting of data in channel i, /y7 



Setting the partial derivatives with respect to a, b, and 
c equal to zero and cancelling common terms gives three 
equations with three unknowns: 



2 4 3 2 

Zw.x. lny. = aZw.x. + bZw.x. + cZw.x. 
li ii ii ii 



3 2 

Zw.x. lny. = aZw.x. + bZw.x. + cZw.x. 
l l i ii ii ii 



Zw.lny. = aZw.x. + bZw.x. + cZw. 

i J i ii ii i 



Applying Cramer's Rule [Ref. 6] to the system determines 
the coefficients: 
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where 

det (A) 



det (A 1 ) 



det (A 2 ) 



det (A^) 
a = det (A) 

det(A 2 ) 
b = det (A) 

det (A^) 
c = det (A) 



= Zw. x. 4 Zw . x. 2 {Zw . - Ew . x . } 

11 i i L i i i J 

3 2 3 

+ Ew.x. {Ew.x. Ew.x. -Ew.x. Zw.} 

i i li i i i l i J 

2 3 2 2 

+ Zw.x. {Zw.x. Zw.x. -(Zw.x. ) } 

li li 1111 

3 

= Zw.x. { Zw . x . Zw . lny . -Zw . Zw.x. lny . } 
i i i i i ^i ill 2 i J 



2 2 2 
+Zw i x i {Zw i x i Zw i x i lny i -Zw i x i Ew i x ± lny i } 

2 2 
+ Ew i x i {Ew i x i Ew i x i lny i -Ew i x i Ew i x i lay^ 



= Ew. x. {Ew. Ew.x. lny . -Ew. x. Ew. lny. } 
xi x xx x xx x -*x 



3 2 2 

+ Ew.x. {Ew.x. Ew . lny . - Ew . Ew . x . lny . } 
X X X* 1 1 X X X X X 



2 2 2 
+ Zw i x i {Ew i x i Ew i x i lny i -Ew i x i Ew^lny^ 
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4 2 

det(A~) = Ew.x. {Ew.x. Ew . lny . -Ew . x . Ew . x . 

■J 11 11 1 1 1111 



3 2 3 

+Ew.x. {Ew.x. Ew.x. Iny. .-Ew.x. Ew . lny . } 
ii 1111 J i ii i J i 



+ Ew^x^^ { Ew^x^^Ew^x^lny^-Ew^x^^Ew^x^^lny^ } 



Exactly the same forms for the coefficients result from 
applying Gauss- Jordan elimination [Ref. 6]. A total of 
eight distinct weighted-data summations form the basic ele- 
ments which combine to give the coefficients. Coefficient 

"a" and a channel number-to-log voltage conversion factor 

2 

deliver the lognormal variance (C^ ) which, with the optical 
wavelength and range, determines the refractive index struc- 
ture constant (C N ) . Coefficients "a", "b” , and the conversion 
factor fix the average log voltage value (In 

C. EXTINCTION 

Extinction is the loss of time-averaged intensity with 
distance by a beam propagating through the atmosphere. No 
extinction measurements had been taken by the optical propa- 
gation group before September, 1975. Assuming an exponential 
dependence and spherical spreading, extinction can be 
represented by an equation of the form: 

I o E o 2 

I = 2 ex P “ R 0 ) ] 



where : 
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